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I discuss recent developments in the QCD cascade formalism. I focus on the 
importance of and uncertainties in higher order corrections to the Modified Leading 
Log approximation for final-state radiation. I also talk about the CCFM and LDC 
evolution equations for initial-state radiation in DIS. 



1 Introduction 

The QCD cascade is an essential part of our understanding of multiparticle 
production in the strong interaction. It is the basis of many tests of per- 
turbative QCD. Furthermore, in order to disentangle more and more subtle 
non-perturbative effects, noting that non-perturbative QCD and the confine- 
ment mechanism is one of the main unsolved problems of high energy physics, 
precise knowledge of the cascade phase is needed. 

Thus, there have been activities aiming at improving cascade 
calculations ,B@|joEll3 and I give a brief status report of these in this talk. 
The content may appear negative and stressing problems, but that is just be- 
cause I want to focus on unsolved questions that hopefully will be the topics 
of progress in the near future. 

2 e + c annihilation 

It is convenient to illustrate the QCD cascade with the logarithmic phase 
space triangle. With transverse momentum and rapidity expressed in terms 
of light-cone momenta, 

k\=k + k^, y=iln(|±), (1) 
we find that an upper constraint \fs/2 on the energy + k-)/2 implies 



Ivl < Mf^) + H 1 + V \ 4kl/s ) « \ []*(■) - Hkl)] ■ (2) 

Thus, the allowed phase space is approximately a triangle in the (y,ln(fc^)) 
plane, cf. Fig [|. This triangle is suitable for QCD discussions, as the gluon 
emission density off a qq dipole behaves like dyd(ln(fcj_)). 
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Figure 1. left: A triangle in the (rapidity, ln(fc^)) plane is a good approximation to the 
phase space given by a maximum energy (dashed hyperbola) and a minimum kj_ (base 
line), right: After a gluon emission, the total base line of the two new dipole triangles is 
larger than ln(s). The additional phase space can be drawn as a double sided fold. 

In the dominant case of significantly different fcj_ , the emission density for 
two gluons off a qq pair is well approximated by one factor representing the 
emission of the first (hardest) gluon g off the qq colour dipole, times another 
factor representing the density for the second (significantly softer) gluon. The 
latter is given by the summed densities off two independently emitting dipoles, 
one for the qg dipole, one for the gq dipole. (There is also a negative term 
suppressed by a relative factor 1/N% which reduces the colour factor in the 
quark directions from N c /2 to Cp, but this will be neglected in this simplified 
discussion.) 

The iterative cascade is based on the observation that this factorization 
generalizes, provided the cascade is strongly ordered in k±. In the parton 
formulation of QCD cascades,!] the emission density off a dipole is split in 
two terms, representing the emission off each individual parton of the dipole, 
respectively. Colour coherence then leads to the well-known angular ordering 
constraint Jj which implies that the parton cascade evolves in angle. 

After the emission of the first gluon g, the qq phase space triangle is not 
large enough to represent the two new dipoles. In the right half of Fig. |l|, 
the additional phase space is drawn as a double sided fold where the two new 
triangles meet. The top of this fold will then have coordinates at 

y s = hn( S ^ l n( fci g) =ln(^M). (3) 

The fold represents the phase space for emissions along the new gluon, with 
rapidities larger than the latest emission, in other words, the angular ordered 
cone around the gluon. For each subsequent emission, a new fold is added, and 
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Figure 2. The parton picture of the cascade evolves in angle (left) and the dipole picture 
evolves in k± (right). 

thus the triangular picture with folds illustrates the iterative QCD cascade. 

Alternative to the parton formulation, the cascade can be interpreted in 
ternitL of dipoles, where each gluon emission splits one dipole into two new 
ones.u The dipole cascade evolves in k±, measured in the rest-frame of the 
emitting dipole, which makes the formalism manifestly Lorenz invariant. As 
illustrated in Fig ||, the parton and dipole formulations of QCD cascades are 
similar, though the evolution parameters differ. Within the accuracy of the 
modified leading log approximation (MLLA), EH they are identical. El 

The MLLA of QCD cascades systematically includes corrections sup- 
pressed by a relative factor Ja^. It is quite possible to go further and system- 
atically include hi gher order corrections embedded in the evolution equations, 
both in the partor£3a and dipoletru formalisms. 

The results from such an exercise must be interpreted with some care. 
Corrections beyond MLLA accuracy depend on the treatment of a very hard 
first gluon, and also "moderately ordered emissions" , by which I mean two 
subsequent gluon emissions with very similar k± . For these, the factorization 
ansatz of the iterative cascade does not hold. Furthermore, one order beyond 
MLLA accuracy enter energy conservation effects. This implies that each 
emitted gluon not only adds phase space to the emitting dipole, but also "eats 
up" phase space in the neighbouring ones. The evolution parameter is then 
promoted, from a convenient book-keeping tool of the cascade, to a physical 
assumption about which gluon indeed is being emitted before others. In a 
sense, we then take the classical cascade formulation of quantum mechanical 
multiparton production more seriously than we are allowed to. 

The study of higher order corrections in cascades is, nevertheless, justified, 
for a very simple reason. Corrections one order beyond MLLA are numeri- 
cally large, for some easily examined observables like inclusive multiplicity 
distributions. 
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Figure 3. Cascade predictions for the ratio of multiplicities in gluon and quark jets, N g /N q 
(left) and the anomalous dimension N g /N g (right). The leading order result (double log 
approximation) is shown with solid line, MLLA results with crosses. Corrections one order 
beyond MLLA are numerically large, and the parton (dashed), dipole (dotted) and gen- 
eralized dipole (dash-dotted) equations differ noticeably. The two dipole alternatives give 
essentially identical predictions for N g /N q , and only the dipole equation results are shown. 



Calculating these corrections, we get an estimate of the important energy 
conservation effects and, by comparing different cascade pictures, we get a 
hint to the uncertainties in these estimates. 

Fig. |^ shows predictions on the ratio of multiplicities in gluon and quark 
jets, N g /N q , and the anomalous dimension N'/N g . Results are obtained with 
three alternative cascade evolution equations. One for the parton picture, and 
two for the dipole picture. The most recent one, called the "generalized dipole 
evolution equation" is presented in detail elsewhere.Q The stronger reduction 
of N q /N g obtained in the dipole equation is favoured by data.E. 



3 Deep Inelastic Scattering 

In Deep Inelastic Scattering (DIS), the magnitude of the phase space triangle 
is determined by the mass W of the hadronic system. The phase space for final 
state emission off the struck quark is given by the photon virtuality scale Q, 
and is represented by a lower corner of the total triangle. A virtual propagator, 
for which ln(fc^) > ln(fc + fc_), can be represented by a line connecting the k+ 
and /c_ values at the k\ height. 

In DIS, high-fc^ emissions, in particular in the target region, are sup- 
pressed. Therefore, the magnitude of the triangle is not as directly linked to 
the multiplicity as in the corresponding case for a qq pair. The picture is 
nevertheless a very useful illustration of DIS events. 

Fig. H illustrates how the cascade for initial-state radiation in DIS starts 
at a cut-off virtuality Qo and an energy fraction xq (selected from an input 
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2ln(Q) ' - iri{l7x) 



Figure 4. The phase space triangle in DIS. The lower left corner represents the phase 
space for final-state emissions off the struck quark. The horizontal line represents a virtual 
propagator with ln(fc^) > ln(fc_|-fc_). 




Figure 5. DGLAP (left) evolution dominates for large ln(Q 2 ) and moderate ln(l/x), while 
non-DGLAP (right) evolution is more important for small ln(Q 2 ) and large ln(l/x). Folds 
related to emitted partons are not shown. 

structure function) and then evolves to the photon interaction point (top 
of photon triangle). For large Q 2 /Qo and moderate 1/x, we have DGLAE 
evolution, where emissions ordered in rapidity are strongly ordered in fc^J13 
For small Q 2 /Qo and large 1/x, this need not be the .case. 

For high enough energies, the BFKL evolution,t3 where emissions are 
strongb: ordered in ln(l/x) but unordered in k±, is expected to dominate. The 
CCFMtl evolution, and the alternative formulation called the "linked dipole 
chain model" (LDC)Q, extrapolates smoothly between the DGLAP and BFKL 
regimes. 

Parton cascades in DIS aim at describing both total cross sections and 
final-state distributions. In the following, I will reduce the final-state distri- 
butions to one observable, which is transverse energy in the target direction, 
called "forward E±" . 

Predictions are obtained from evolution equations convoluted with input 
structure functions at some cut-off scale Qq. Fig. |6| shows what happens when 
Qo is raised. The DGLAP evolution predicts no significant forward E±, and 
this simple prediction is of course independent of Qo- A change in Qo is 
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Figure 6. As Fig. q but with a higher value of the cut-off scale Qoi which implies that 
perturbative calculations are not trusted in the gray region. In the non-DGLAP case 
(right), this gives events whose contribution to forward E± is not easily included in the 
formalism. 

compensated by a change in the input parton distributions, leaving the total 
cross section unchanged. 

Non-DGLAP evolution predicts more forward E±, but a quantitative pre- 
diction may be difficult to make comfortably insensitive to the cut-off scale 
Qq. With a rise in Qq, some previously allowed propagators of the evolu- 
tion may drop below the cut-off, which in principle implies that we no longer 
trust their behaviour to be determined by perturbative physics. Thus, it is 
not clear how to properly treat emission chains like the one above the gray 
region on the right hand side of Fig. ^, and thus how to properly include their 
contribution to forward E± . 

More investigations are needed to answer whether this contribution is nu- 
merically important or not. It may be that the discussion above is a problem 
in principle only, and not in practice, but it is interesting to note that Monte 
Carlo investigations of the CCFM and LDC schemes find forward E± to be 
highly sensitive to subleading corrections in the evolution, and that these cor- 
rections move predictions too far below data, while leading order calculations 



4 Summary 

As a natural consequence of the success of the parton cascade formalism, 
higher order corrections are being examined. These are found to be 

• theoretically uncertain, 

• numerically important. 

This should of course be interpreted like this: QCD in Multiparticle pro- 
duction is an interesting field of physics, suitable for cooperation between 
experimentalists and theorists! 
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